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Abstract 

It has been proposed that the energy evolution of QCD amplitudes 
in the high-energy regime falls in the universality class of reaction- 
diffusion processes. We review the arguments for this correspondence, 
and we explain how it enables one to compute analytically asymptotic 
features of QCD amplitudes. 

The high-energy regime of QCD has been intensively studied in deep-inelastic e-p scat- 
tering at HERA, in heavy-ion collisions at RHIC, and will be probed in proton and heavy-ion 
scattering at LHC. It is rich of interesting theoretical structures: Links have been found or con- 
jectured with conformal field theory, string theory, and more recently, with a class of models 
known in statistical physics. 

When hadrons scatter at very high energy, the color fields that are generated at the in- 
teraction point have a large strength. Perturbative methods qualify as soon as there is a large 
transverse momentum scale in the event: This property enables one, under certain conditions, 
to derive QCD evolution equations, in the form of partial differential equations (which can also 
be stochastic in particular formulations such as the dipole model [1]). On the other hand, strong 
fields cause the parton densities to saturate, which makes this evolution nonlinear. 

Similar looking stochastic nonlinear partial differential equations also appear in problems 
of apparently different physical origin, such as reaction-diffusion, or population evolution. The 
goal of this short review is to explain that these similarities are not casual, and that once under- 
stood, they can help the derivation of new results for QCD cross sections. We refer the reader to 
the original papers [2, 3] for the details, and to Ref. [4] for a more extensive review. 

In the following, we will consider the scattering of two hadrons, and we will aim at com- 
puting their cross section at very high energies. Their relative rapidity is denoted by Y. Since our 
discussion will rely on resummed perturbative QCD, we think of these hadrons as being small 
objects, such as color dipoles found for example as fluctuations of highly virtual photons. We 
will always be discussing a definite region of impact parameter. 

1 High energy QCD and reaction-diffusion 

Cross sections are measured by counting the number of events that are registered in a detector 
within a given interval of time. Each single event results from an interaction between the scatter- 
ing hadrons realized as definite quantum states, that is, as particular Fock states. Let us go to the 
frame in which one hadron is almost at rest, while the other one carries most of the kinetic energy, 
and thus develops a highly-occupied Fock state. As long as saturation effects are negligible (i.e. 
far from the unitarity limit in which the hadrons appear black to each other), the probability of 
interaction is proportional to the number of partons in the fast hadron whose transverse momenta 



k match the typical momentum scale of the slow hadron. Let us imagine that k is tunable (it 
is the case when the slow object is a virtual photon), and that one could actually measure the 
interaction probability of the slow hadron with a particular Fock state of the fast hadron. (In 
practice, it would require the replication of quantum Fock states, which is impossible). We call 
T{k) this interaction probability (more precisely, it is the forward elastic scattering amplitude at 
a fixed impact parameter). T is an unphysical quantity, but it will be important to understand it 
theoretically. The physical amplitude A(Y, k) is just the average of T{k) over all possible Fock 
state realizations at the considered rapidity Y: A = (T), that is, A is the average over all events 
that may occur, appropriately weighted by their probabilities at a given rapidity. 

Standard quantum field theory calculations, based on the evaluation of Feynman diagrams, 
would directly lead to the expression of A. However, it turns out that such calculations are ex- 
tremely hard. Instead, understanding first the main analytical features of the scattering probability 
T off a single typical Fock state of the fast hadron and then averaging over events is a simpler 
approach, that has been successful in leading to analytical expressions for the asymptotics of the 
scattering cross sections. 

What is precisely known about T is its evolution with rapidity (or energy), at least in the 
regime in which T< 1, that is, away from the unitarity limit. When one increases infinitesimally 
the rapidity of a hadron that has say n partons in its current Fock state, there is a transition rate to 
a n + 1, n + 2...-parton Fock state that is computable in perturbative QCD. It may be extracted 
from the BFKL equation. A direct formulation is the color dipole model [1], in which this 
transition probability is explicitely computed. When T <C 1, T is a linear function of the number 
of partons. Roughly, it reads T(k) ~ a 2 s n(k), where n(k) is the number of gluons that have a 
transverse momentum of the order of k. 

This transition to higher Fock states may be captured by a linear stochastic equation, of 
the form 

da Y T = X (-d hlk 2)T + a s VTu, (1) 

where a = a s N c /ir. v is an appropriate stochastic variable that has zero mean, and variations 
of order unity when aY is increased by one. \{— d\ Ilk 2) is the usual BFKL kernel. It describes 
the branching diffusion of partons, at least in the regime of very high energy in which we are 
interested in. This means that, when acted on T, it roughly behaves like a diffusion term df nk2 T 
supplemented by a growth term T (with appropriate coefficients). The noise term is a conse- 
quence of discreteness: It implements the fact that we are considering the evolution of one single 
Fock state, that contains a definite (discrete) number of partons. 

On the other hand, when T becomes of the order of 1, saturation effects have to enter in 
order to tame the growth of the number of partons, for unitarity to be preserved. From the work 
of Balitsky and Kovchegov (BK), we know that in the mean-field limit in which the noise can 
be neglected (that is, when A = T; This is realized when one of the interacting objects is a very 
large nucleus), the evolution equation for T reads 

d 5tY T = X (-d hlk 2)T-T 2 . (2) 



Hence we shall propose that the full evolution be described by the following stochastic equation: 

d a YT = x(-d hlk 2)T-T 2 + a s VTv. (3) 




Fig. 1: Sketch of the gluonic content (the gluons are represented by disks) of the hadron in the transverse plane. In a 
little boost, each gluon can split, or nonlinear effects like recombination can take place. At the time of the interaction, 
the system is probed by a dipole of size l/k, which is sensitive to the number of gluons of similar size at a given 
impact parameter. 

This equation is in the universality class of the stochastic version of the Fisher and Kolmogorov, 
Petrovsky, Piscounov (F-KPP) equation. (The latter would in fact be obtained by replacing 
x(— d\ n k 2 )T by di nk 2T + T and \fT by y/T{l — T); for a review and references, see [5].) 
A schematic picture of the evolution is presented in Fig. Q] 

So far, no one has succeeded in formulating rigorously scattering in QCD in the form of 
a stochastic equation such as Q, in particular for the way how saturation occurs is not yet fully 
understood. Is it gluon recombination, as was advocated in the early papers on saturation [6], or 
some more subtle process? So the best one can do at this stage is to set the noise term in such 
a way that the evolution of the hadron Fock state by gluon splittings is exactly reproduced away 
from the unitarity limit, where the nonlinear term may be neglected. A practical implementa- 
tion of this process would be, for example, Salam's Monte-Carlo code of the dipole model [7] 
modified by the addition of a suitable saturation condition which makes sure that T (written in 
coordinate space) keeps always less than 1 and thus that unitarity is preserved. In this procedure, 
the expression for the noise is unambiguously fixed (it results from the splitting probability of 
the dipoles), and the BFKL limit is exactly taken into account. Such a procedure was suggested 
in Ref. [8], but has not been implemented so far for its technical awkwardness. Some other 
paths were followed: One may alternatively take v to be a Gaussian white noisd3 [10]. This 
simple choice enables one to apply the Ito stochastic calculus, and to draw a link with equations 
established within QCD such as the B-JIMWLK equations (For a review and references, see 
Ref. [11], and A. Shoshi's talk at this conference [12].) One may also think of the whole process 
as a reaction-diffusion process, as we will implicitely do in the next section. 

All this may look quite arbitrary: We have merely merged two known limits into a single 
equation, without much further justification. How can one be sure that A obtained from aver- 
aging realizations of Eq. ([3]> looks like the solution of a genuine QCD equation? Although it 

'However, in this case, T should not exactly be the amplitude, but rather a kind of "dual amplitude" - see for 
example Ref. [9]. 



might sound strange a priori, the solution to the "arbitrary" equation that we have written down 
is very likely to contain the exact asymptotics of QCD. This fact is actually related to the univer- 
sality of the solutions to such evolution equations. The statement is the following: For a large 
class of processes, i.e. for a number of stochastic functions v and for a variety of forms of the 
nonlinearities, the asymptotics of the statistics of T (that is, the physical observables A = (T), 
(T 2 )...) for small a s and large aY are identical. This is not a theorem, but a conjecture based on 
a general understanding of how noisy traveling waves propagate. (The propagation mechanism 
is described in the next section.) The whole point is that the details of the evolution equation do 
not matter for extracting the asymptotics of the QCD amplitudes. 

Let us describe a typical process whose evolution is in the universality class described 
by the F-KPP equation: reaction-diffusion. This process involves particles on a lattice indexed 
by some variable x, that evolve by a set of rules of the following form: As time is increased, 
each particle has a probability either to jump to a nearby site, or to become two particles, or 
to recombine with another particle on the same lattice site. The balance between creation and 
recombination of particles determines the equilibrium number of particles on each site N. After a 
large evolution time, the number of particles on a given site oscillates about N (with an amplitude 
of the order of the typical statistical fluctuations y/N). It is the number of particles per site 
normalized to N that obeys an equation in the universality class of the F-KPP equation. 

At this point, we may establish a simple dictionary between reaction-diffusion and QCD. 
Time is the evolution variable, so is rapidity: t <-> aY. The variable in which diffusion takes 
place is x <-> In k 2 . The equilibrium number of particles is N <-> 1/a 2 . (In QCD, it is fixed by 
the unitarity conditiorJ^T < 1.) 

2 Statistical methods and application to QCD 

In a first step, we ignore the stochastic term, that is, we address the simple BK equation (O, in 
order to gain intuition on the form of the solutions. A given localized initial condition (T ~ a 2 s in 
a region of order 1 around some initial scale In k 2 : This would be the physical initial condition) 
will spread and grow under the action of the kernel \{~ d\nk 2 )> which, as we wrote before, 
amounts to a branching diffusion. But as soon as T becomes of the order of 1, the nonlinear 
term enters to compensate the growth, making T saturate. Then further evolution necessarily has 
the form of two symmetric traveling waves, since the system can only escape to the right and to 
the left. Let us focus on the rightmoving wave, that travels towards larger values of In A: 2 . This 
wave front is represented schematically in Fig. [2] It turns out that the shape of this wave in its 
large-ln k 2 tail is exponential, with a slope that is completely fixed by the linear kernel: 

T ~ e" 70 ln fc2 , 70 being determined by = X ' (70) • (4) 

7o 

2 This condition actually holds in coordinate space (when T is a function of transverse sizes). In momentum space, 
the growth of T with energy is also tamed as soon as the point T = 1 has been crossed, although T(k) can take 
arbitrarily large values. This does not change the conclusions that we shall draw later: The only important feature of 
the evolution is that T changes behavior in the saturation region. One can see how it goes precisely in QCD e.g. in 
the numerical simulations presented in Ref. [13]. 
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Fig. 2: Deterministic F-KPP front and its evolution with time. The arrows show where branching diffusion takes place 
to drive the motion towards larger values of x. 

Since the wave front keeps its shape, it makes sense to characterize its motion by a single velocity 
Voo- The latter is also completely determined by the kernel x- It reads 

Voo = *M. (5) 
7o 

In QCD, the position X t of the wave is called the saturation scale Q S (Y). It characterizes the 
momentum below which nonlinear saturation effects (gluon recombination, multiple scattering...) 
become important. The velocity defined above is simply the derivative of h\Q 2 (Y) with 
respect to Y. (Recall that the x- variable is Ink 2 .) 

Now that we have understood the deterministic limit, we may try to put back the noise. We 
do not know how to attack the problem directly. Instead, we shall stick to a physical approach, 
and view the evolution equation (0) as describing a reaction-diffusion process. In this framework, 
we recall that the origin of the noise was the discreteness of the number of particles on each 
site. Discreteness means in particular that the number of particles n(k) cannot be a fraction of 
an integer. Consequently, coming back to the simple-minded relationship^ T(k) ~ a 2 n(k), it 
means that T is either are larger than a 2 . Brunet and Derrida [15] proposed to replace the full 
stochastic equation by a deterministic one that takes into account this basic effect of discreteness, 
which can easily be done by not allowing any growth when T < a 2 . (It amounts to cutting 
off the tail of the front; to do this in practice, one may for example replace \ by a modified 
kernel obtained by subtracting its growth term in the region in which T < a 2 . Note that there 
is no unique prescription.) The solution to this modified equation is again a traveling wave, 
that exhibits the same overall exponential decay as given by Eq. dU) (except for an uninteresting 
additional pref actor). Its velocity now read^l 

T/ x(7o) t 2 7ox"(7o) 
7o 2ln 2 (l/a 2 ) 



3 Again, this is not literally true: T(k) is actually continuous, but the tails (below T = a 2 ) are decaying exponen- 
tially with a characteristic length of one unit in the variable In k 2 . This is steep enough for all our arguments to apply 
as if T(k) itself were discrete. 

4 This result had already been obtained by Mueller and Shoshi [14]. Actually, the understanding of high-energy 
scattering as a peculiar reaction-diffusion process emerged from a reinterpretation of their work, in the light of the 
Brunet-Derrida analytical treatment of traveling waves with a cutoff [15]. 




Fig. 3: Traveling waves solution of the F-KPP equation with a cutoff that simulates discreteness. In the left sketch, 
no stochasticity is taken into account while in the right one, particles may randomly be sent ahead of the deterministic 
front. Their further time/rapidity evolution is also represented. 

It is thus less than the velocity of the front in the limit of an infinite number of particles (obtained 
by letting a s go to 0), which had to be expected: Indeed, taking into account discreteness amounts 
to removing some "matter" from the front, which logically slows down its motion. 

However, a deterministic solution can only reproduce approximately the realizations of a 
stochastic evolution. We can incorporate stochasticity back into the picture [16] by noting once 
again that the noise is only important in the forward tail of the front, where the number of particles 
is low on the average. From numerical simulations of simple reaction-diffusion models, we 
observed the following behavior: Most of the time, the motion of the front is almost deterministic, 
with a velocity given by the solution to the cutoff deterministic equation. From time to time, 
rarely, a large fluctuation causes a transitory acceleration of the front. This fluctuation consists in 
one or a few particles being sent far ahead of the deterministic tip of the front, which then evolve 
into a new front that later gets absorbed by the deterministic front. This behavior is represented 
in Fig. [3] We conjectured a probability distribution for these fluctuations, as well as the effect 
that they have on the position of the front after relaxation. 

With these elements, we were able to deduce the full statistics of the saturation scale, that 
is to say not only the mean position (or velocity) of the front, 

(lnQ 2 s ) X (7o) ?r 2 7 o X "(7o) 2 2 31njn(l/a|) 
V = -v = n1 27i / 9x + n ^oX (7o) — ; 37- , 9X , (7) 



but also all its cumulants: 



7o ln d (l/a2) 



when n > 2. 



Now we recall that the physical amplitude is obtained by averaging T over all possible 
realizations. Given that the fall off of the large-ln k 2 tail of each single event is exponential, it is 
not difficult to get the scaling of the scattering amplitude with the help of Eq. ©: 

ln a (l/"§) 




Fig. 4: Time/rapidity evolution of a noisy traveling wave. The noise is essentially concentrated at the tip of the front, 
where the occupation numbers are low. After some time (rapidity), the wave has moved to the right (3 realizations 
are shown with thin lines), roughly keeping its shape. However, stochasticity manifests itself macroscopically by 
inducing a dispersion in the positions of the fronts between different realizations. Since the physical amplitude A is 
the average of all realizations, its very shape is influenced by the noise. {A is represented by the thick line.) 

where (In Q 2 S (Y)) is given by Eq. f7]). This is the main analytical result for QCD that comes out of 
the statistical approach. Note that other results can be extracted on the statistics of the branchings 
of the gluons in the course of the evolution, but we cannot see a possible phenomenological 
application. 

The emerging overall picture of front propagation is shown in Fig. @] 
3 Prospects 

Clearly, the statistical interpretation of scattering processes has proved useful since it has led to 
both a new understanding and new asymptotical results for high energy QCD. Of course, it relies 
on a few conjectures that will eventually have to be proved in a more formal way, but we feel that 
we have so far provided robust physical arguments. 

It has to be acknowledged that our new analytical results are not relevant for phenomenol- 
ogy yet, since they make sense for hx{l/a s ) 3> 1 only, which requires values of a s so small 
that, of course, they are far beyond the experimentally attainable range. A number of authors 
have however taken seriously the extrapolation of these results to realistic values of a s and have 
produced predictions, see e.g. Ref. [17]. On the other hand, numerics could give results valid 
for a s < 0.1 (optimistically), that is, not far from the phenomenological domain. (This point is 
discussed in Ref. [8]). 

At this point, we have been able to extract properties that QCD shares with simple sta- 
tistical models. We could claim that this was a correct procedure because asymptotic properties 
of the solutions do not depend on the details of how saturation occurs. So in some sense, we 
have done the "easy" part of the work. However, to go closer to phenomenology, one would 



need to understand more deeply the details of saturation, which probably constrain the form of 
the noise v in Eq. ((3]). Investigations of some possible models have been conducted, sometimes 
leading to pecularities in the interpretation, such as negative transition rates [18]. Building a 
complete picture, valid beyond asymptotics, remains a challenging open question, for which a 
further breakthrough may be needed. 

Finally, our approach to the propagation of noisy traveling waves is not based on a field 
theory formulation, but is an event-by-event analysis of the shape of realizations, using methods 
more familiar to statistical physicists than to particle physicists. Being able to recover results 
such as Eqs. ([71), ®, © within field theory, starting e.g. from an effective Lagrangian whose 
building blocks are Reggeon fields, would be a very interesting achievement. Some progress has 
been made recently, see e.g. Ref. [19]. 
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